IR Phonon signatures of multiferroicity in TbMn 2 Os 
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The infrared (IR) active phonons in multiferroic TbM^Os are studied as a function of tempera- 
ture. Most of the symmetry allowed IR modes polarized along the a and b crystal axes are observed 
and the behavior of several b polarized phonons is correlated with the magnetic and ferroelectric 
transitions. A high frequency b polarized phonon only Raman allowed in the paraelectric phase be- 
comes IR active at the ferroelectric transition. The IR strength of this mode is proportional to the 
square of the ferroelectric order parameter and gives a sensitive measure of the symmetry lowering 
lattice distortions in the ferroelectric phase. 

PACS numbers: 78.30.-j 63.20.Ls 75.50.Ee 



I. INTRODUCTION 

Materials that exhibit simultaneously (anti)- 
ferromagnetic, (anti)-ferroelectric or (anti)-ferroelastic 
degrees of freedom, the so called multiferroics, have 
become of increasing interest because they offer great 
possibilities for applications in multifunctional devices. 
These applications are specially based on the power 
to control the magnetic state by electric fields and 
the ferroelectric state by magnetic fields, this type of 
materials being the magnetoelectrics, a subclass of the 
multiferroicsi2ii4. Moreover, the understanding of the 
basic mechanism that allows this behavior is important 
in the development of the devices that make use of these 
characteristics. 

The interplay of these degrees of freedom has been 
demonstrated recently in several multiferroic materials 
such as i?Mn 2 05 where R is Eu, Gd, Tb, Dy, Ho and 
Y 5 i 6 i 7 i 8 i 9 as well as in TbMnO^ and Ni 3 V 2 0g±I. Fun- 
damental questions as to the nature and characteristics 
of the coupling of the magnetic and ferroelectric orders 
remain unanswered. Phenomenological approaches" 
based on the Landau theory of phase transitions have 
been used to relate the symmetry of the spin order with 
the appearance of ferroelectricity. The basic conclusion 
of these proposals is that the non-collinearity of the spin 
system is crucial in the development of ferroelectric or- 
der. More microscopic studies^i^ have also pointed to 
the importance of non-collinearity of the spin system. 

The antiferromagnetic manganese oxide TbMn 2 05 (or- 
thorhombic space group Pbam # 55, Z = 4) is a multi- 
ferroic with non-collinear magnetic orde»i& that shows a 
strong magnetoelectric coupling effect: an applied mag- 
netic field along the a axis changes the sign of the electri- 
cal polarization present along the b axis. The dielectric 
constant e along the b axis has anomalies associated with 
the distinct phase transitions at low temperatures (see 
figure 1 in Hur, et al 6 ): at the Neel temperature Tn ~ 
42 K no anomalies are present in e. There is a paraelec- 
tric to ferroelectric phase transition at Tc ~ 38 K evi- 



dent by a peak in the dielectric constant. The magnetic 
order then locks in to a commensurate structure (CM) 
with wave vector (1/2,0,1/4). At T w 24 K the mag- 
netic order transforms into an incommensurate structure 
(ICM) with a step-like feature in e; this anomaly is also 
accompanied by hysteresis^. We note that the dielectric 
constants along the a and c axis show no significant ef- 
fects at these phase transitions. 

In EuMn 2 Os Polyakov, et aU£ found displacements of 
the Mn 3+ ion along the a axis at the ferroelectric transi- 
tion. They suggested that this behavior leads to a change 
in symmetry from the space group Pbam to the non- 
centrosymmetric group P2iam 26). In a analogous 
work on the compound YM^Ojr 1 - 8 , Kagomiya, et al pro- 
posed similar displacements at the ferroelectric transi- 
tion. The displacements reported are very small (« 0.007 
A) , which hints to a different origin of the ferroelectric- 
ity when compared to the typical proper ferroelectrics. 
Other structural worksi^^i 9 ^ have not reported any sig- 
nature of atomic displacements at the ferroelectric phase 
transition in this family of compounds. Moreover, Mi- 
hailova, et al2i reported a Raman study in H0M112O5 
and TbMn 2 5 and Garci'a-Flores, et alS 2 . in BiMn 2 5 , 
EuMn 2 Os and DyMn 2 Os as a function of temperature 
and found no evidence of anomalous behavior of the Ra- 
man active phonons at the ferroelectric transition tem- 
perature. Nevertheless, the existence of a macroscopic 
dipole moment is evidence of the lack of inversion sym- 
metry in the FE phase, and motivates the study of the 
dynamics of the lattice to look for further information 
about the structural changes. 

In this report we present a study of the temperature 
dependent infrared (IR) phonon spectra of the multifer- 
roic TbMn 2 05. The most interesting feature we find is 
the appearance of an IR inactive phonon activated at the 
ferroelectric transition with light polarization parallel to 
the static ferroelectric polarization Po (l? w ||Po||i>). This 
indicates that one IR forbidden mode (Raman or silent) 
in the paraelectric phase acquires an electric-dipole mo- 
ment due to the static displacement associated with the 
ferroelectricity. We identify this phonon with a Raman 
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FIG. 1: (Color online). Experimental (solid line) and fit 
(dashed line) at T = 7 K phonon spectra in TbM^Os. 
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FIG. 2: (Color online). Optical conductivity of the newly 
activated phonon (shifted plots) . Arrows indicate the position 
of the resonance frequency. 



Mn-0 stretching mode, which accounts for its sensitiv- 
ity to the static polarization possibly induced by the Mn 
spin systeniiiSi2£. 



II. EXPERIMENTAL RESULTS 

Single crystals of TbMx^Os were grown using 
B 2 3 /PbO/PbF 2 flux in a Pt crucible. The flux was 
held at 1,280 °C for 15 hours and slowly cooled down 
to 950 °C at a rate of 1 °C per hour. Crystals grew 
in the form of black platelets as well as cubes with 
a typical size of 10 mm 3 with a working diameter of 
5 mm. The crystals were characterized and oriented 
using x-ray diffraction at room temperature. Normal 
incidence reflection spectra were taken with a Bomem 
Fourier Transform Spectrometer DA3.02. Light was po- 
larized along the a and b crystal axes in the frequency 
range 8-800 cm -1 (w 1-100 meV) and in the tempera- 
ture range 7 to 300 K. Sample was kept in vacuum in 
a continuous He flow cryostat with optical access win- 
dows. The factor group analysis based on structural 
data by Alonso, et a,& of the paraelectric phase gives 
the following IR active vibrational modes at the T point: 
T IR = 8B lu (E\\z) + UB 2u {E\\y) + UB 3u (E\\x) identical 
to a previous report^i. We complemented this analysis 
with a shell model calculation of the phonon dispersion. 

Figure ^ shows the spectra at T = 7 K. Twelve of 
the 14 IR active phonons in the a polarization were re- 
liably observed, whereas all 14 phonons polarized along 
b are present in the spectrum. The reflectivity spectra 
was fitted in a least squares procedure using the sum of 
Lorentzian form of the model dielectric function e, given 



by: 

N Si 
e(uj) = £oo + V— i (1) 

where is the dielectric constant at high frequency 
and the phonon parameters ujq, S and 7 are the phonon 
frequency, the spectral weight and the linewidth respec- 
tively, we also define Si = Ae^o;^ where Ae is the con- 
tribution of the phonon to the static dielectric function; 
these parameters are extracted as functions of temper- 



TABLE I: Oscillator parameters at Ti = 7 K and T 2 = 45 K 
in TbMn 2 5 . a, b are the crystal axes. e^ b = 5.31, 6.82. 
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FIG. 3: (Color online). Comparison of the behavior of the 
spectral weight and frequency of the new phonon to the static 
polarization Po. Polarization and frequency show warming 
curves and S shows the cooling curve. 



ature and are displayed in table |U The result of the 
fitting is also shown in figure ^ and it can be seen that it 
is almost indistinguishable from the data indicating the 
weakness of higher order phonon processes. 

We obtained the optical conductivity using the 
Kramers-Kronig transform of the reflectivity spectrum. 
Figure shows the optical conductivity around 700 cm -1 
for several temperatures with _E w ||Po. We observe that a 
feature not present at 45 K appears in the low T phases. 
The temperature dependence of the spectral weight and 
frequency of this feature are plotted in figure where 
we see the spectral weight starting to appear at 38 K. 
S for this phonon was obtained by directly integrating 
the optical conductivity between 695 and 710 cm -1 . The 
spectral weight of this phonon increases and its frequency 
shifts, both continuously as we lower the temperature. 
Around 24 K both abruptly change and show hysteresis 
around this point evident by the difference in the cool- 
ing and warming curves. The behavior of this feature is 
correlated with the second order FE transition at 38 K 
and the first order CM — > ICM transition at 24 K in this 
compound. The static polarization P plotted in figure|3| 
was obtained by measuring the temperature dependence 
of the pyroelectric current on a similar sample to the one 
used in the optical measurements. 

On very general grounds we can relate the appearance 
and behavior of this phonon to these underlying phase 
transitions. Since the lattice distortions Su associated 
with these phase transitions are very small we can ex- 
pand the spectral weight and frequency shifts in pow- 
ers of Su. The quadratic term is the first non-zero term 
in this expansion that can describe the spectral weight 
change or frequency shift. Similarly the order parame- 
ters associated with the new phases are proportional to 
Su so that Pq oc Su. As a result we expect that the 



TABLE II: Irreducible representation splitting at the FE 
phase transition. S. A. = Spectral Activity (R = Raman 
active, IR = Infrared Active). 



# 55 S. A. 


# 55 irreps 
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#26 S. A." 


R 


A a 


Ai 


IR (y) & R 


R 


B ig 


B 2 


IR (x) & R 


R 


B 2g 


A 2 


R 


R 


B Zg 


B x 


IR (z) & R 


Silent 


A u 


A 2 


R 


IR(z) 


B\ u 


Si 


IR (z) & R 


m(y) 


B 2 u 


Ax 


IR (y) & R 


IR(x) 


Bz u 


B 2 


IR (x) & R 



"Note that in this column x, y, z correspond to the high temper- 
ature system of coordinates a, b, c and differ from what is found 
in the character table. 



spectral weight behavior and frequency shift should be 
S, Auj oc (Su) 2 oc Po 2 - This is the observed behavior as 
can be seen in figure where we have plotted Po 2 with 
the phonon data. At low temperatures (T < 10 K) where 
the Tb moment orders, the phonon data deviates from 
Po 2 suggesting that it is the Mn and oxygen ion displace- 
ments that dominate the dynamics of this high frequency 
phonon. The possible scenarios for the appearance of a 
new phonon are: (1) zone folding of the phonon disper- 
sion (since the magnetic order corresponds to a lock in 
ICM -> CM transition with k = (1/2,0,1/4)), and (2) 
activation of IR-inactive phonons at this transition due 
to the loss of inversion symmetry. Our shell model cal- 
culation shows that the dispersion of the high frequency 
phonon is negative so that no zone-folded mode can give 
the high frequency of this phonon. We therefore conclude 
that this phonon is a previously IR inactive phonon that 
acquires electric dipole moment at the FE transition. 

In a ferroelectric phase transition, where inversion 
symmetry is lost, symmetry considerations dictate that 
phonons that were not IR active in the paraelectric phase 
can become IR active in the FE phase. This is the case 
in TbMn2 05, where the low T phase has mixed IR and 
Raman phonons. In this low T phase the phonons of 
the high T symmetry group split as shown in table ITT? 4 . 
This splitting was obtained by considering what sym- 
metry operation is maintained in both phases and then 
assign spectral activity according to the experimental ob- 
servations. The assumption was made as well that the 
low T space group is as proposed by Polyakov, et alii 
and Kagomiya, et ati&. Consistency with the laboratory 
frame (x — > a, y — > b, z — > c) is applied as well. From the 
reports by Mihailova, et ali2i and Garcia-Flores, et al^2 
we learn that a A g mode at frequency rj 700 cm -1 ex- 
ists in all the RM.YI2O5 materials whose Raman phonons 
have been reported. We note as well that these reports 
do not resolve any IR phonons becoming Raman active at 
the FE phase transition. We conclude that this high fre- 
quency A g Raman phonon is the mode we observe that 
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FIG. 4: Temperature dependence of the b axis Tb phonon fre- 
quency. Inset shows the full temperature dependence. Typ- 
ical dispersion of several measurements is indicated as error 
bars. 



FIG. 5: (Color online). Temperature dependence of the spec- 
tral weight of 2 oxygen dominated phonons. Inset shows the 
temperature dependence of the spectral weight of one of the 
Tb 3+ crystal field transitions. 



acquires IR activity in the FE phase. Further experi- 
mental confirmation in the Raman spectrum is needed 
to make a definitive identification. 

Only a few other phonons show correlations with the 
low temperature phase transitions. The phonons polar- 
ized along the a axis do not show any significant anoma- 
lies in this temperature range. This is consistent with 
the featureless behavior of the dielectric function along 
this axis. On the other hand, the behavior of some of 
the phonons with dipole moment along b is non-trivial. 
The low frequency phonon with frequency w 96 cm , 
identified primarily with movement of the Tb ions, has 
a temperature dependence that correlates with the low 
temperature CM — > ICM magnetic transition. In figure 
0|the frequency of this phonon is plotted versus tempera- 
ture and we observe an increase in the frequency around 
24 K. This effect is thought to be a manifestation of the 
coupling of this phonon to a magnon as is discussed by 
Katsura, et al2£. Further experimental results on this 
effect will be presented elsewhere^. 

Surprisingly several phonons show interesting temper- 
ature dependence for T above Tjv- The inset in figure 0] 
shows the full temperature dependence of the frequency 
of the b axis phonon. The anomalous softening in the 
temperature range of 150 K to 50 K demonstrates addi- 
tional effects in the dynamics of the lattice. In figure^the 
behavior of the spectral weight of two oxygen phonons 
polarized in the a and b axes, with frequencies of 704 
and 689 cm" 1 respectively, seems complementary: the a 
phonon gains spectral weight while the b phonon looses it 
as the temperature is lowered. This effect is present in the 
full temperature range from 300 K to 7 K, while for the 
rest of the phonons the spectral weight is only changed 
significantly around the various transition temperatures. 
The fact that these modes gain or lose so much spectral 



weight (10 and 6-fold respectively) in a large tempera- 
ture interval also demonstrates some higher energy scale 
in this system. These effects (the anomalous softening 
and the dramatic changes in spectral weight) are not un- 
derstood at present. However, one interesting possibility 
is that they are of magnetic character. Recent high tem- 
perature susceptibility measurements 22 in BiMx^Os has 
shown evidence for spin frustrated behavior from devia- 
tions from Curie law with a Weiss temperature of s» 250 
K. Dielectric anomalies in B1M112O5 around this temper- 
ature have been reported 2 ^ as well. 

Finally the inset of figure [S] shows the temperature de- 
pendence of the intensity of a feature observed at 120 
cm" 1 . This is identified as a crystal field level of the Tb 
3+ ion 28 . This transition has electric dipole character as 
is seen from the form of the reflectivity curve (see fig. 
as well as the fact that the spectral weight (see ta- 
ble P) is comparable to the IR active phonons (magnetic 
dipole transitions are usually much weaker than electric 
dipole transitions). This conclusion is supported as well 
by the shell model calculation that shows the 3 lowest 
phonon excitations being the Tb-dominated phonon (at 
tzs 100 cm" 1 ) and then a doublet (at rj 170 cm" 1 ). Fur- 
thermore, the observed temperature dependence of the 
intensity is common for the f-level crystal field transi- 
tions in the rare earth ionsSS. Several other absorption 
features are observed in transmission at lower frequencies 
and will be reported elsewhere 2 ^. 



III. CONCLUSIONS 

We have measured the IR phonon spectra in TbMn 2 05 
along the a and b axes and have observed most of the 
symmetry allowed modes. The majority of the phonons 
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do not show significant correlations to the FE and AFM 
phase transitions of the system. However several phonons 
exhibit interesting correlations to the ferroelectricity of 
this material. We have found a signature of the loss of 
inversion symmetry in the FE phase by the appearance 
of a IR phonon below T c that was only Raman active 
in the paraelectric phase. The strength of this mode is 
proportional to the square of the FE order parameter 
and gives a sensitive measure of the symmetry lowering 
lattice distortions in the ferroelectric phase. In addition, 
the lowest frequency phonon (along b) displays hardening 
in the CM — > ICM transition possibly due to the coupling 
with the spin system^; 2 modes (along a and 6) have 
dramatic changes in their spectral weight over a wide 



temperature range possibly because of frustration effects 
in the spin system. We have also identified an electric 
dipole active crystal field transition of the Tb 3+ ion in 
the phonon frequency range. 
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